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ABSTRACT 



After the initiation of the explosion of core-collapse supernovae, neutrinos emitted from the nascent neutron star drive a supersonic 
baryonic outflow. This neutrino-driven wind interacts with the more slowly moving, earlier supernova ejecta forming a wind termi- 
nation shock (or reverse shock), which changes the local wind conditions and their evolution. Important nucleosynthesis processes 
(alpha-process, charged-particle reactions, r-process, and vp-process) occur or might occur in this environment. The nucleosynthesis 
depends on the long-time evolution of density, temperature, and expansion velocity. Here we present two-dimensional hydrodynami- 
cal simulations with an approximate description of neutrino-transport effects, which for the first time follow the post-bounce accretion, 
onset of the explosion, wind formation, and the wind expansion through the collision with the preceding supernova ejecta. Our results 
demonstrate that the anisotropic ejecta distribution has a great impact on the position of the reverse shock, the wind profile, and the 
long-time evolution. This suggests that hydrodynamic instabilities after core bounce and the consequential asymmetries may have 
important effects on the nucleosynthesis-relevant conditions in the neutrino-heated baryonic mass flow from proto-neutron stars. 

Key words, supernovae: general — neutrinos — nuclear reactions, nucleosynthesis, abundances — hydrodynamics 



1. Introduction 

Supernova outflows are an important astrophysical site for sev- 
eral nucleosynthesis processes. A variety of isotopes seem to be 
exclusively produced there and contribute to the metal enrich- 
ment of the interstellar medium from which old halo stars and 
later our Solar System have formed. Therefore, the fingerprints 
of supernova nucleosynthesis are searched for in the old halo 
stars and are atte mpted to be extracte d from the Solar System 
abundances (e.g., ISneden et all I2008I) . Studying the long-time 
evolution of core-collapse supernovae is a challenging prob- 
lem b ecause the explo sion mechanism is not yet fully under- 
stood ( Jan ka et all I2007 ) and because of the high computational 
demands for simulating such a dynamical environment. CPU 
time requirements become even more extreme when multidi- 
mensional simulations combined with accurate neutrino trans- 
port are supposed to follow the supernova ejecta for many sec- 
onds. In this paper, we take advantage o f a computationally ef- 
ficient neutrino transport approximation (iScheck et all I2006I) to 
present the first results of two-dimensional simulations of core- 
collapse supernovae that track the outflow evolution for up to 
three seconds. 

Simultaneously with the onset of the core-collapse super- 
nova explosion, the proto-neutron star forms and cools by emit- 
ting neutrinos (iPons et all Il999l) . The latter deposit energy in 
the layers near the proto-neutron star surface and thus drive 
a baryonic mass outflow. Because of ongoing neutrino-energy 
input, the entropy of this wind grows with distance from the 
neutron star, and the pressure gradient exceeds the gravita- 
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tional force of the compact remnant. Therefore the neutrino- 
heated matter is accelerated quickly and even reaches super- 
sonic v elocities (Duncan et al], Il986l) . After first promising re- 
sults of IWoosley et al.l ([1994), neutrino-driven winds have been 
investigated intensely as a site where hea vy elements could be 
produced via rapid neutron capture (see lArnould et all 120071 
for a review of r-process sites and physics). However, the 
condit ions for successful r-processing found by IWoosley et al.l 
(1994) could neither be confirmed by other sim ulations at that 
time (Witti et all Il994t iTakahashi et all Jl994). nor by mor e 
recent ones ( e.g., | Arcones et all 12007k iFischer et all 120101) . 
iRoberts et alJ (12010) ' explain the numerical problems that arti- 
ficially produced s uitable conditi o ns for the production of heavy 
n-rich elemen t s (iHoffman et all 1 19971) in the simulations of 
IWoosley et al.l (1 19941) . Because of the lack of the appropriate 
conditions in more recent neutrino-driven wind simulations, ef- 
forts continue to find possible missing physical ingredients. In 
lArcones et al.l (120071) . hereafter Paper I, we studied in detail the 
evolution of the neutrino-driven wind and its interaction with 
the earlier supernova ejecta for different stellar progenitors and 
neutrino luminosities. This interaction results in a wind termi- 
nation shock (or reverse shock), which changes the evolution 
of nucleosynthesis-relevant conditions: density, temperature, ex- 
pansion velocity (see Paper I). Howe ver, integrated nuc l eosyn- 
thesis based on those simulations ( Arcones & MontesL l2010h 
show that no heavy r-process elements can be produced. 

First systematic studies of neutrino-driven winds were based 
on solving stationary wind equations (see e.g.JOian & Woosleyl 
Il996l:l0tsuki et alll200d ; lThompson et alll200ll) . These steadv- 
state models could not consistently account for the reverse 
shock, which is a hydrodynamical feature found in several 
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(e g.,lJanka& Mulled. 1995; Janka & Miiller, 
al., 1995; BurasetaL, 2006; Arcones et al, 



1996; 



2007 



simulations i 
Burr ows et < 

iFischer etall 120101) . The impact of the reverse shock on the r- 
process has been investigated by means of parametric, steady- 
state models of subsonic "breeze" solu tions combined with an 
outer boundary at c onstant pressure (iSumiyoshi et all [2000; 
iTerasawa et all 120021) or supersonic winds with fixed as ymptotic 
temperature (e.g., IWanajo et all 120021: IWanajol 120071) . or two- 
stage outflow models whe re a faster wind phase is followed by 
a slow er expansion phase dKuroda et all 120081; iPanov & Ja nka, 
20091). There are also recent studies (IWanajo et all 12010: 
Roberts et alll2010|) of the effect of the reverse s hock on the vp - 
process (Frohlich et al ll2006l:IPruet et all 120061: 1 WanaioLl2006l) . 
A discussion of the implications of the reverse shock on the nu- 
cleosynthesis will be given in Sect, ffl 

In Paper I, we have analyzed the behaviour of the reverse 
shock based on spherically symmetric hydrodynamic simula- 
tions. We found by analytic means that the position of the re- 
verse shock depends on wind velocity and density, but also on 
the pressure of the supernova ejecta with which the wind col- 
lides. This pressure is strongly related to the explosion energy 
and progenitor structure. Therefore, the evolution of a mass ele- 
ment that has crossed the reverse shock is more complicated than 
a simple boundary at constant pressure or temperature. Here we 
show that multidimensional effects lead to an anisotropic mass 
and density distribution of the initial supernova ejecta, which has 
a big impact on the position of the reverse shock. The location 
of the latter becomes angle dependent. In our two-dimensional 
(2D) simulations, the neutrino-driven wind boundary is therefore 
strongly deformed to a non-spherical shape. As a consequence, 
the properties of the shocked wind medium are strongly depen- 
dent on the direction. The results presented here will help to im- 
prove the simple extrapolations for the ejecta evolution that are 
used in nucleosynthesis studies. They will allow to constrain the 
parameter space of possible wind histories. 

The paper is organized as follows. Numerical method and 
computed models are described in Sect. O Our results of two 
simulations for a chosen stellar progenitor are presented in 
Sect. where we also compare 2D to ID results (Sect. 13.21) and 
analyze the impact of varying the progenitor (Sect. 13.3b . A de- 
tailed discussion of the assumptions, approximations, and future 
improvements is included in Sect. 13.41 Finally, the possible nu- 
cleosynthesis implications of our results are addressed in Sect.H 
and we summarize our findings in Sect. 



2. Two-dimensional simulations 

Observations indicate that core-collapse supernova are highly 
anisotropic. Therefore, multidimensional simulations are more 
realistic than spherically symmetric ones. However, they are 
computationally much more expensive. Using the same hydro 
code as in Paper I, we have performed hydrodynamical simu- 
lations of the neutrino-driven wind phase in two dimensions. 
We follow the evolution of the supernova ejecta for two sec- 
onds starting at a few milliseconds after bounce. As in Paper I, 
our simulations are done with a Newtonian hydrodynamics code, 
which incl udes general relativ ity corrections in the gravitational 
potential (iMarek et all 120061) . It is combined with an efficient 
neutrino transport approximation (IScheck et all [2006). In order 
to improve the performance of our simulations, the central high- 
density region of the neutron star is excised and its behaviour is 
prescribed by an inner boundary condition. This allows us also to 
vary the contraction and final radius of the neutrons star, which 
are determined by the uncertain high-density equation of state. 



Neutrino luminosities at the boundary are chosen such that the 
explosion energy is around 1.5 B. A detailed description of the 
numerical method, in itial models, a nd b oundary treatment can 
be fo und in Paper I, IScheck etall (120061) . and iKifonidis et al.1 
d2003l) . 

Here we discuss six different models: two (T10-ll-rl and 
T10-12-rl) are based on a 10.2 M Q star (data provided by 
A. Heger, priv. comm.), three (T15- 12-rl, T15-ll-rl, and T15-1 1- 
rO) on a 15 M© progenitor (si 5 s7b2. 1 Woosley & Weaver] (fl99 5)), 
and on e (T25-15-r4) on a 25 M© progenitor (s25a28. LHeger et all 
( 2001)). All progenitors were followed through core collapse by 
A. Marek and are mapped to our 2D code typically 10 ms af- 
ter bounce. In all models the subsequent contraction of the neu- 
tro n star and its coolin g behaviour are described as explained 
by IScheck et al.l (120061) and defined by parameters Rf for the 
final boundary radius and to for the contraction timescale (see 
Table []])• The neutrino luminosity imposed at the lower grid 
boundary is + L % = 52.5, 38.6, and 30.3 B/s for the mod- 
els including "11", "12", and "15" in the name, respectively. For 
the model names we follow the same convection as in Paper I. In 
Table [T] we summarize the values of the input parameters. The 
initial configuration of the progenitor model is spherically sym- 
metric. Since the code conserves this symmetry, it is necessary 
to add nonradial perturbations to some quantity to trigger the 
growth of nonradial hydrodynamic instabilities in regions where 
the conditi ons allow for the d evelopment of such instabilities. 
Following IScheck et al.l (120061) we apply seed perturbations of 
the velocity field in our simulations, for which we choose cell- 
to-cell random variations with a typical amplitude of 0.1 %. The 
resolution of our two-dimensional simulations is around 900 ra- 
dial grid points and 1 80 angular bins. The number of radial zones 
is increased depending on the requirements in the outer layers of 
the proto-neutron star, where the density gradient is very steep 
and neutrinos decouple from matter. 

3. Results 

In this section we describe the explosion and ejecta evolution 
for models T10-ll-rl and T10-12-rl, give an analytic explana- 
tion of the reverse shock behaviour in Sect. 13.11 and compare to 
spherically symmetric models (Sect. 13.2b and to the more mas- 
sive progenitors (Sect. l3.3l) . The model parameters and overview 
of the results are presented in Table [T] The proto-neutron star 
contraction is characterized by its time scale to and the final 
boundary radius Rf. The boundary luminosity (L* + Lf) is con- 
stant during the first second and later decays as r 3/2 . The end of 
the simulation is denoted by the time £ e nd given in seconds after 
bounce. Mb ar is the baryonic mass of the neutron star. The neu- 
tron star radius R ns is defined as the location where the density 
is 10 11 gem -3 . A£ tot is the total energy radiated in neutrinos of 
all flavors (measured in bethe [B] = 10 51 erg), L Ve and L Ve are the 
luminosities of electron neutrinos and antineutrinos at 500 km, 
(£ Ve ) and (£ Ve ) are the corresponding mean energies, E ex p is the 
explosion energy, t Qxp is the post-bounce time when the explo- 
sion sets in (defined as the moment when the energy of expand- 
ing postshock matter exceeds 10 49 erg). 

The initial distribution is spherically symmetric except for 
small random seed perturbations. As neutrinos deposit energy 
behind the shock, a negative entropy gradient establishes. The 
region between neutron star and shock thus becomes Ledoux- 
unstable and convective overturn appears. Neutrino-heated, high 
entropy matter streams upwards, while downflows transport low 
entropy matter from the shock to the neutron star. The down- 
flows and rising bubbles evolve quickly on time scales of the 
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Model 


Contraction 




Progenitor Mass 


^end 


M bar 


A£ tot 


Rns 






fee) 


fee) 




^exp 




(Rf, to) 


[B/s] 


[M ] 


[s] 


[M ] 


[100B] 


[km] 


[B/s] 


[B/s] 


[MeV] 


[MeV] 


[B] 


[s] 


T10-ll-rl 


9km; 0.1s 


52.5 


10 


2.8 


1.261 


1.305 


14.82 


22.97 


24.63 


20.51 


22.10 


1.457 


0.153 


T10-12-rl 


9km; 0.1s 


38.6 


10 


2.2 


1.280 


1.146 


13.44 


21.76 


22.49 


21.36 


22.91 


0.938 


0.170 


T15-12-rl 


9km; 0.1s 


38.6 


15 


1.5 


1.421 


1.228 


12.76 


22.60 


23.23 


22.27 


23.75 


1.405 


0.184 


T15-ll-r0 


8km; 0.1s 


52.5 


15 


2.0 


1.393 


1.460 


12.79 


25.53 


26.45 


22.68 


24.04 


1.364 


0.156 


T15-ll-rl 


9km; 0.1s 


52.5 


15 


1.0 


1.388 


1.461 


13.27 


27.66 


28.22 


22.55 


23.87 


1.341 


0.162 


T25-15-r4 


10.5 km; 0.1s 


30.3 


25 


1.6 


1.869 


2.233 


13.41 


36.46 


39.92 


24.31 


25.51 


3.674 


0.197 



order of tens of milliseconds. The mass distribution below the 
shock becomes highly anisotropic and is dominated by low 
spherical harmonics modes. In this phase continuous neutrino 
heating aided by convection leads to an explosion. The evolu- 
tion after the onset of the explosion is shown in Fig. [T] by the 
entropy distribution at different times after bounce. 

After the launch of the explosion a proto-neutron star forms 
at the center and cools and deleptonizes by emitting neutrinos. 
During this phase, neutrino-heated matter expands away from 
the proto-neutron star surface, in what is known as neutrino- 
driven wind and collides with the previous, more slowly moving 
ejecta. The interaction of the supersonic wind with the super- 
nova ejecta results in a wind termination shock (reverse shock). 
At this discontinuity, kinetic energy is transformed into internal 
energy, which produces an increase of the entropy (see panels of 
Fig. [T] at t > 500 ms) and the wind is decelerated. At late times 
the changes become slower, shock and ejecta expand quasi- self- 
similarly. 

The neutrino-driven wind is spherically symmetric because 
it depends only on the neutrino emission of the proto-neutron 
star, which is isotropic. However, the distribution of the early 
ejecta is highly anisotropic. This produces a deformation of the 
reverse shock and of the shocked wind. The reverse shock radius 
and the properties of the shocked matter become angle depen- 
dent. An example of a two-dimensional feature is the downflow 
present at 6 « n/2 in model T10-ll-rl, which corresponds to the 
low entropy region visible in the panels for t = 500 - 2000 ms 
of Fig. [T] This feature leads to a big deformation of the reverse 
shock, in contrast to model T10-12-rl, where the reverse shock 
stays almost spherically symmetric as no strong, long-lasting 
downflows have formed. These variations in the evolution and 
structure of the two models are not an immediate consequence of 
the different boundary luminosities. In Table Q] one can see that 
many parameters of both models are very similar. The anisotropy 
of the ejecta depends on chaotic effects that are trig gered by the 
initial random perturbations (for more details see IScheck et all 
120061) . 

High-density, low-entropy regions in the ejecta, which are 
the remainders of former downflows, act like obstacles prevent- 
ing faster wind expansion in those directions. In the analytic dis- 
cussion of Paper I, we found that the reverse shock radius (R YS ) 
depends on the pressure of the more slowly moving early ejecta: 



R r „ oc 



The aspherical matter distribution in the layer between re- 
verse shock and forward shock leads also to an angle depen- 
dence of the entropy jump, which depends on the reverse- shock 
position and wind v elocity (approximately as s TS oc ^/R^ul^ 5 , 
I Arcones et al.l (120071) ). Figure [3 shows the radius of the reverse 
shock (bottom panel) and the entropy, pressure, radial velocity, 
and temperature of the decelerated wind matter just after the re- 
verse shock as functions of the azimuthal angle for model T10- 
11-rl at t = 1.5 s after bounce. Pressure and reverse shock radius 
are related as roughly given by Eq. (Q]). Moreover, the pressure 
determines also the angle between radial direction and reverse 
shock (Sect. 13.11) . The temperature, which is one of the most 
relevant quantities for nucleosynthesis, depends mainly on the 
position of the reverse shock. When the reverse shock is at a 
smaller radius the shocked matter reaches higher temperatures. 
The link between these quantities and the geometrical structure 
of the reverse shock can be explained by analytic means as we 
show in the next section. 



3.1. Analytic discussion 

We use the Rankine-Hugoniot conditions (see paper I for further 
discussion) to understand the behaviour of the shocked matter 
and its angular dependence. The neutrino-driven wind expands 
in the radial direction with velocity w w , hits the slow-moving, 
anisotropic ejecta, and a deformed reverse shock results. As 
shown schematically in Fig. [3j there is an angle between re- 
verse shock and radial direction. For spherical explosions this 
angle is always n/2, since the reverse shock can only be perpen- 
dicular to the radius vector. In an oblique shock the velocity can 
be decomposed into two components: tangential (u\\) and perpen- 
dicular (w ± ) to the shock. The tangential com ponent of the veloc- 
ity is continuous through an oblique shock ( Landau & Lifshitz, 
1 1959b . Following the notation introduced in Fig.O this implies 
that: 

u\\ = u w cos = w rs?r cos + u Y ^e sin (2) 



with 



-*rs,r 



Urs,e = w rs smx , 



(3) 
(4) 



(1) 



being the radial and lateral components of the velocity of the 
shocked matter (w rs ), respectively. 

The perpendicular component of the velocity, u ± , is changed 
according to the Rankine-Hugoniot conditions. The mass con- 
servation condition, including the angle dependence, is: 



The mass outflow rate (M w ) and the velocity (u w ) of the wind are 
the same at all angles because the wind is spherically symmetric. 
Therefore, the variation of the reverse shock radius with angle is 
caused by the pressure variations of the anisotropic ejecta CP rs )- 



p w w w sin0 = p rs w rs , r sin0 + p rs w rs ^cos0. (5) 
The momentum continuity condition for an oblique shock is: 

P w + p w u 2 ± = P rs + Prs4,± > ( 6 ) 
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Fig. 1. Entropy distribution in models T10-ll-rl (left column) and T10-12-rl (right column) for different times after bounce as 
indicated in every panel. The figures are plotted such that the polar axis is oriented horizontally with "south" (0 = n) on the left and 
"north" (0 = 0) on the right. The thin grey line marks the shock radius. In the panel for t = 1500 ms of model T10-ll-rl, the radial 
lines mark the angular directions at # = 25 degrees (green line) and 100 degrees (red line), along which radial profiles are shown in 
Fig. El 
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Fig. 2. Latitudinal variation of the reverse shock radius (bottom) 
and of the entropy, pressure, radial velocity, and temperature of 
the wind matter just after passing the reverse shock in model 
T10-ll-rl at t = 1.5 s after bounce. 



where only the perpendicular components of the velocities enter. 
We can thus write a relation between the pressure and the angle 




Fig. 3. Schematic representation of the velocities in a fluid going 
through an oblique shock (red line). Only the velocity compo- 
nent perpendicular to the shock, u±, is changed when a mass 
element crosses the shock; the tangential component, u\\, is con- 
served. Therefore, the direction of the flow is changed at the 
shock. 



AP = p w ^sin 2 ( 



i-I 



(7) 



where AP = P rs - P w 
(note that usually P w 
(Fig.O, and/3 is: 

/3 = P -H = 



is the pressure jump at the reverse shock 
«: Prs), u w is the radial wind velocity 



u± 



u w sin 



w rs?r sin + w rSj cos 



(8) 



The variation of the pressure jump (and thus of the entropy and 
temperature increase by the reverse shock) with angle can be 
deduced from Eq. ©. Since J3 > 1 usually, this equation implies 
that when goes to n/ 2, i.e., the reverse shock is perpendicu- 
lar to wind velocity, AP and the entropy jump are larger (see 
middle panel in Fig.©. As there is (roughly) a pressure balance 
between P rs and the pressure of the slow-moving early ejecta, 
which is higher in the downflows, tends to be about n/2 in 
regions where downflows are present. In Fig. [H the upper panel 
shows that the density reaches highest values where the down- 
flow is located (6 « 100 degrees). Consistently, the pressure in 
this region is also large as shown in the middle panel. This leads 
to straight sections in the reverse shock shape and to the occur- 
rence of kinks between locations of downflows. 

The reverse shock in model T10-ll-rl thus exhibits several 
kinks (Fig. [4]) due to the anisotropic pressure distribution of the 
ejecta. Figure \5\ shows, in a simplified way, the effects of these 
kinks. An oblique reverse shock is less effective in decelerating 
the flow (Eq. ([8])) because the tangential component of the veloc- 
ity is conserved through the shock. This leads to higher veloc- 
ities outside the non- spherical parts of the reverse shock (u YS in 
Fig. 0) compared to spherical regions (u' rs in Fig.0. This can be 
seen in the velocity field marked by arrows in the middle panel 
of Fig. |4l The velocity of the ejecta is not directed radially after 
the wind has passed an oblique shock. This produces collimated 
high- velocity outflows starting at the kinks (Fig.|4l bottom panel 
and Fig.Q]). Moreover, the kinks of the reverse shock are asso- 
ciated with minimum values of the pressure and entropy for the 
shocked material as visible in Fig.[2l 
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Fig. 4. Distribution of density, pressure, and absolute value of 
the velocity in the wind and shocked matter for model T10-ll-rl 
at t - 1.5 s after bounce. In the middle panel the radial velocity 
field is indicated by arrows. 




Fig. 5. Schematic representation of the velocities in a fluid going 
through the reverse shock (red line) when a kink leads to the 
collimation of the outflowing matter. 



3.2. Two-dimensional versus one-dimensional simulations 

Convection enhances neutrino heating, leading to earlier and 
more energetic explosions in 2D than in ID, for the same inner- 
boundary parameters. Earlier explosions imply that the neutron 
star accretes matter during less time, thus its mass is slightly 
smaller in 2D (see Tabled . These differences alter the wind and 
reverse shock properties, making it difficult to contrast exactly 
one- and two-dimensional simulations. In this section we com- 
pare the one-dimensional model M10-ll-rl (see Paper I) to the 
two-dimensional model with the same inner boundary parame- 
ters: T10-ll-rl. First, we examine the differences between the 
radial profiles of both models at a given time and later the evo- 
lution of relevant quantities. 

Figure shows the radial profiles of different quantities at 
1.5 s after bounce. The one-dimensional model M10-ll-rl is dis- 
played by black lines. For the two-dimensional model, T 10-11 - 
rl, the radial profiles differ at all azimuthal angles. Therefore, 
we chose two angles (see radial lines in panel for t = 1500 ms 
of Fig. [T]): = 25 degrees (green line) and 6 = 100 degrees (red 
line), where R YS adopts extreme values. As in Paper I, the radial 
profiles show an increase of the entropy in the region where neu- 
trinos deposit energy (r < 20 km) and a constant entropy value 
in the wind. Matter close to the neutron star absorbs neutrinos 
and moves outwards. The fast drop of the density and pressure 
allows the expansion to become supersonic. The wind velocity 
increases first approximately as u oc r and then approaches an 
asymptotic value (w w ). The interaction of the fast wind with the 
slow-moving earlier ejecta results in a sudden drop of the ve- 
locity and a jump of the pressure, density, temperature, and en- 
tropy to higher values at the reverse shock. While the wind is still 
very similar in the one- and two-dimensional simulations, there 
are significant differences in the conditions of the slow-moving 
ejecta and therefore in the reverse shock position. The variations 
in the wind are just a consequence of different neutron star prop- 
erties due to the evolution towards explosion, which is not the 
same in one and two dimensions, as explained before. However, 
the two profiles of model T10-ll-rl are identical in the wind 
phase indicating that this region is spherically symmetric. This 
is expected because the neutrino emission is isotropic and the 
neutron star stays spherical without rotation. The differences in 
the profiles appear at the position of the reverse shock, which in 
2D depends on the angle and is also different from the ID case. 
Notice that for = 100 degrees the velocity of the slow-moving 
ejecta is very small because of the presence of relics of a strong 
downflow. 

Since the main differences between ID and 2D originate 
from the reverse shock radius and are linked to the properties 
of the supernova ejecta, we show in Fig. [7] the evolution of the 
quantities of matter that has passed the reverse shock. The evolu- 
tion in 2D is presented for fixed azimuthal angles = 25 degrees 
(green lines) and 100 degrees (red lines), corresponding to the 
lines of same colors in Fig. [T] In 2D the reverse shock radius 
changes significantly between different angles with correspond- 
ing variations of the temperature. This has an impact on nucle- 
osynthesis as we will discuss later in Sect. [4] The behaviour of 
the reverse shock radius follows Eq. (Q]). The entropy in the two- 
dimensional simulation is always lower, even when the reverse 
shock radius is larger than in the one-dimensional case. This is 
explained by the fact that the reverse shock for the green ray is 
oblique to the radial direction and the entropy jump is thus re- 
duced (see Sect. 13.11) . Moreover, in 2D the neutron-star mass is 
slightly smaller, which le ads to lower entropies in the wind itself 
dOian& WoosleyLll996h . 
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Note that the wind phase can be studied just by one- 
dimensional simulations because (without rotation) it stays 
spherically symme t ric. Even simple steady - state wind models 
lOtsuki et al] (120001) : lThompson et al.l (l200ll) are sufficient to de- 
termine the wind properties, e.g. for discussing the nucleosyn- 
thesis in the wind. However, the interaction of the wind with the 
supernova ejecta is a hydrodynamical problem that requires su- 
pernova explosion simulations. Moreover, we have shown here 
that one-dimensional models are not enough to account for all 
the possibilities of the long-time evolution of the ejecta. 

3.3. Progenitor dependence 

The reverse shock depends on the pressure of the more slowly 
moving earlier ejecta, which is different for different progenitor 
stars, as shown in Paper I. For similar explosion energies, more 
massive progenitors have slower ejecta and therefore the ejecta 
shell possesses a higher pressure so that the reverse shock stays 
at a smaller radius. We have seen in the previous section that an 
anisotropic distribution of the pressure in 2D has a big impact 
on the reverse shock position. It is thus interesting to study the 
combined effect of the two ingredients: progenitor structure and 
anisotropic ejecta. 

In Fig. [8] the entropy distributions are shown for the rest of 
the models in Table [TJ representing the explosions of 15 M Q 
and 25 M© progenitors. Note that the time is different in every 
panel because the panels correspond to the moments £ e nd when 
our simulations were stopped. We can compare models with the 
same parameters for the proto-neutron star evolution, but dif- 
ferent progenitors, e.g., T 10-11 -rl and T15-ll-rl. The evolution 
of the ejecta in these two models is considerably different, al- 
though the boundary conditions are the same. This is even more 
visible in the case of models T10-12-rl and T15-12-rl. The wind 
termination shock is almost spherically symmetric in the 10 M© 
star, while it is highly asymmetric due to the presence of long- 
lasting downflows in the 15 M© progenitor. The third model for 
the 15 M© case (Fig. [8j upper panel) also develops a relatively 
spherical reverse shock. Although there is no unambiguous re- 
lation between the evolution of the ejecta and the shape of the 
reverse shock on the one hand and the boundary conditions or 
the progenitor structure on the other, more massive progenitor 
stars like T25-15-r4 favor long-lasting and more slowly expand- 
ing downdrafts because of their denser structure and higher mass 
infall rates. The same trend can be observed in the cases of lower 
boundary luminosities (e.g. in T15-12-rl) and less energetic ex- 
plosions. For the explosion models of the 10 M© star shown in 
Fig.[U the expansion is relatively fast because of the more dilute 
silicon and oxygen shells, and therefore all downflows are blown 
away during the first half a second, even when lower luminosi- 
ties are assumed in the simulations. 

However, the ejecta distribution does not only depend on the 
progenitor and on variations of the boundary conditions that in- 
fluence the strength of the explosion. Initial random perturba- 
tions must be imposed on the progenitor (or post-bounce) model 
to trigger the growth of instabilities and anisotropics . These de- 
velop in a stochastic and chaotic way and thus can lead to signif- 
icant variations of the ejecta morphology, i.e., of the number and 
position of downflows and the direction and size of high-entropy 
bubbles, all of which have an influence on the deformation of the 
forward and reverse shocks and thus on associated variations of 
the nucleosynthesis-relevant conditions in the early supernova 
ejecta and in the neutrino wind. 

In spite of that, for roughly similar explosion energy the pro- 
genitor structure has a systematic influence -at least in a con- 



strained sense- on the evolution of the supernova shock and of 
the wind-termination shock: As seen from the comparison of 
model T10-ll-rl (Fig. ffl left) with models T15-ll-r0, T15-11- 
rl, and T15-12-rl (Fig. [8]), all of which possess similar explo- 
sion energy, both shocks expand significantly more slowly in the 
case of the more massive progenitor. This trend was already ob- 
tained for ID models in Paper I and is independent of the ex- 
plosion asymmetries that have developed as a consequence of 
the (stochastic and chaotic) growth of nonradial instabilities in 
the shell enclosed by the forward shock and the reverse shock. 
Such a systematic dependence of the ejecta evolution on the 
progenitor structure is therefore independent of the dimension 
and for this reason also of the particular initial seed perturbation 
by which the ejecta asymmetry was initiated to grow in our 2D 
models. The reason for this progenitor-dependent variation is the 
denser structure of the shells around the iron core in more mas- 
sive progenitor stars as discussed in the second paragraph of this 
section and in more detail in Paper I. 

The multi-dimensional situation, however, allows for an even 
more extreme effect that is not possible in the ID case: In pro- 
genitors significantly more massive than our 10 M model (or, 
alternatively, in the case of very weak explosions), long-lasting 
accretion by the neutron star can take place and the explosion 
can become extremely asymmetric such that the neutrino-driven 
outflow can be confined to narrow angular wedges (see models 
T15-12-rl and T25-15-r4 in Fig. [8]). While neutron- star accre- 
tion is still going on, the outflow is then prevented to become 
supersonic in some directions so that a termination shock can 
be absent at these angles. It is therefore clear that basic trends 
defined by progenitor properties (like those mentioned above) 
are superimposed in the multi-dimensional situation by poten- 
tially very large nonradial explosion asymmetries, which cannot 
be predicted merely on grounds of the progenitor structure and 
the value of the explosion energy, but can have a high relevance 
for the nucleosynthetic output of the supernova. 

3.4. Discussion of assumptions, approximations, and future 
improvements 

In the present work our goal was to identify differences between 
the nucleosynthesis-relevant conditions of the neutrino-driven 
wind in two dimensions compared to spherical symmetry. We, 
however, did not intend to cover the topic of multi-dimensional 
effects in the winds in its full breadth, but concentrated in partic- 
ular on the interaction of the wind with the more slowly moving, 
dense ejecta shell behind the supernova shock. This interaction 
has obviously stronger and more important consequences than a 
variety of other nonradial perturbations that can make the wind 
outflow anisotropic at typically much smaller level. 

The wind termination (or reverse) shock, which separates 
the supersonic wind region from the subsonic ejecta shell, ex- 
hibits strong non-spherical deformations because of the large- 
scale inhomogeneities imprinted on the structure of the ejecta 
by the nonradial hydrodynamic instabilities at the beginning of 
the explosion. When the reverse shock is located at a sufficiently 
small radius, it alters the wind conditions at such a high tem- 
perature that the nucleosynthesis in the outflow can be strongly 
affected. Different from the one-dimensional case, where the 
reverse- shock radius is simply a function of the wind properties 
(mass-loss rate and expansion velocity) and explosion properties 
(energy and shock radius) as discussed in Paper I, the location of 
the reverse shock in the multi-dimensional case also depends on 
the inhomogeneities that are present in the ejecta shell. Dense 
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relics of long-lasting accretion funnels, for example, propagate 
more slowly and prevent the reverse shock from rapid expansion. 

In order to gain insight into such effects, we followed 
two-dimensional models for sufficiently long evolution periods 
with the emphasis on maintaining high spatial resolution from 
the steeping density gradient at the proto-neutron star surface 
through the neutrino-wind region, reverse shock, and inhomo- 
geneous ejecta shell, out to the expanding supernova shock. We 
employed a number of assumptions and approximations which 
we will critically discuss in the following. 

We emphasize that in the present work we want to highlight 
basic and general aspects of anisotropic wind-ejecta interaction. 
However, we neither intend to make predictions for particular 
progenitors nor do we want to give a statistical overview of pos- 
sible outcomes. Both are not possible in view of our incomplete 
knowledge of the exact explosion mechanism and the use of two- 
dimensional models. 

3.4.1. Random initial seed perturbations 

In order to initiate the growth of nonradial flows in hydrody- 
namically unstable regions, we perturbed our spherical initial 
models with random seeds (Sect. [2]) . Random perturbations are a 
simple, flexible, and partly tested (see below) procedure to intro- 
duce small- amplitude deviations from spherical symmetry in any 
available progenitor or post-bounce model. Of course, it would 
be interesting to start with a more realistic assumption based 
on predictions from self-consistent stellar evolution modeling. 
However, multi-dimensional simulations were performed only 
for special evoluti on stages of selected progenitor stars (e.g., 
I Young et aill2Q05[ and references therein). Therefore, there are 
no multi-dimensional progenitor data at the onset of the gravita- 
tional instability for a range of stellar masses and the fluctuations 
of quantities in the core at this stage are not well constrained. 

Moreover, tests suggest that the exact properties of the ini- 
tial perturbation pattern are not overly important except in the 
improbable and special case that a low-mode, global asymme- 
try is present in the progenitor core. The influence of different 
random seeds — different amplitudes, different random realiza- 
tions, and also initial perturbations of differe nt quantities — ha s 
been explored (at least to some extent) by iBuras et aT (2006); 
IScheck etall (120061) . One might suspect that a high-mode, low- 
amplitude pattern of the initial perturbations may have an influ- 
ence on the fastest growing mode(s) in the hydrodynamically un- 
stable regions after core bounce as well as on the linear growth 
timescales of these most unstable mode(s). However, it seems 
unlikely that any memory of such a high-mode seed pattern is 
retained once the hydrodynamical instabilities have reached the 
fully nonlinear regime because the involved growth is character- 
ized by highly chaotic and stochastic interactions of structures 
on different scales. 

Bur as et all (|2006) tested the amplification of random pertur- 
bations on the level of 1-2% imposed on the density structure of 
the progenitor core before collapse. They found that these ini- 
tial seeds increase by factors of a few in the supersonically in- 
falling layers during the collapse (and are damped in the s ubson - 
ically falling regions) as predicted by lLai & Goldreichl (2000). 
However, the changes of the model evolution after core bounce 
turned out to be insignificant compared to the standard case in 
which they imposed perturbations on the level of 0.1% just af- 
ter core bounce and thus only shortly before the hydrodynam- 
ically unstabl e conditions develo p in the supernova core (see 
Appendix E in lBuras et aD (120061) ). 



IScheck etaD (120061) investigated the effects of different ran- 
dom seeds imposed on models after core bounce, perturbing dif- 
ferent quantities — radial velocity or density — with an ampli- 
tude of 0.1 percent. They found that different random realiza- 
tions (keeping the progenitor star and boundary conditions the 
same) lead to differences of the convective patterns and thus of 
the detailed structures of the anisotropics in the dense shell be- 
hind the supernova shock. But while the number of long-lasting 
downflows, their locations, the shock deformation and the an- 
gular distribution and volume of buoyant high-entropy bubbles 
varied, little sensitivity to the initial seeds was seen for the 
global properties of the supernova explosion, i.e., the explosion 
timescale and the average shock radius and explosion energy as 
functions of time. Of course, the transition from a long-lasting, 
asymmetric accretion phase to the spherical neutrino wind will 
depend on the size and persistance of the asymmetric accretion 
funnels that are formed when vigorous, non-radial instabilities 
stir the ejecta during the early explosion stage. 

For the present paper we picked a set of models that provides 
a kind of representative overview of different possible morpholo- 
gies (see Figs. [J and [8]) that were also seen in the larger set of 
calculations bv lScheck et al.1 (120061) . 

3.4.2. Anisotropic neutrino emission and flavor conversions 

We describe the neutrino emission from the dense core of 
the proto-neutron star by choosing time-dependent (spherically 
symmetric) luminosities for neutrinos and antineutrinos of all 
flavors at the inner grid boundary, which is placed significantly 
below the neutrinosphere. 

Naturally, the neutrino flux evolution imposed there has a 
strong influence on all stages of the evolution of neutrino-driven 
supernovae: on the shock revival and explosion, on the poten- 
tially long-lasting post-explosion accretion, and in particular 
also on the power and duration of the neutrino-driven wind. 
The dependence of the dynamical and thermodynamical prop- 
erties of the neutrino wind (mass-loss rate, entropy, and expan- 
sion timescale) on the luminosities and mean energies of the ra- 
diated neutri nos was in detail discuss ed in pr evious works, in 
parti cular by lOian & Woosleyl (Il996l) and by iThompson et al.l 
(l200ll) . while the dependence of the wind termination shock on 
the neutrino- wind properties was the subject of Paper I. Making 
assumptions about the neutrino fluxes at the inner grid boundary 
implies that we can only explore possible dynamical scenarios 
in a parameterized way, but we are unable to determine what 
exactly happens in a particular progenitor star. 

Similarly, the neutrino-flavor evolution necessarily plays a 
crucial role for setting the neutron-to-proton ratio in the wind. 
Thus it decides about the nucleosynthesis processes that can take 
place in this environment. While we think that it is not possible 
to make reliable predictions of the electron fraction Y e in the 
wind on the basis of our transport approximations, we are con- 
fident that the set of models presented here can still shed light 
on fundamental aspects of anisotropies that can occur during the 
dynamical evolution driven by the neutrino heating between the 
neutron- star surface and the supernova shock. 

Asymmetric neutrino emission, e.g. associated with non- 
spherical accretion or convective activity inside the nascent neu- 
tron star, might also play a role for creating anisotropies in the 
wind ejecta. The most extreme case of asymmetric accretion are 
long-lasting downflows on one side of the proto-neutron star, 
while mass loss in the wind is going on in other directions. Such 
a situation can — at least temporarily — occur in our simula- 
tions. Because of the violent impact of the accretion flows on 
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the neutron star surface, associated with rapid motions of the 
tips of the accretion funnels, strong sonic waves are launched 
and create non-radi al per turbations in the neutrino-heated out- 
flows (IScheck et al.Ll2Q06h . 

These hydrodynamical perturbations of the winds are usu- 
ally much bigger than any asymmetries caused by anisotropic 
neutrino emission from the accretion downflows. This is so be- 
cause the supersonically infalling matter in the downflows has 
only little time to experience neutrino losses: It reaches suffi- 
ciently high temperatures and densities only shortly before it im- 
pacts on the neutron-star surface. After deceleration, it quickly 
spreads sideways around the surface of the neutron star and set- 
tles into the hot accretion layer surrounding the accretor. While 
the anisotropic emission from the accretion flows contributes to 
the neutrino luminosity only on the level of some percent, most 
of this anisotropic emission occurs in regions where the optical 
depth for neutrinos is much lower than unity. The reabsorption 
of these neutrinos is therefore very improbable and is unable to 
cause any significant asymmetries in the baryonic wind outflow 
(typically such asymmetries remain well below one percent). 
Instead, the wind is driven by neutrino interactions in regions 
much closer to the neutrino sphere, composed of accreted matter 
that has experienced efficient mixing upon being dispersed in a 
turbulent vortex layer. In addition, the strong surface gravity of 
the neutron star prevents asymmetries in the accretion layer so 
that this layer acts as the base of a nearly isotropic wind. 

A convective region below the neutrino sphere, which is only 
partly included in our computational volume (its inner part is 
omitted with the excised neutron-star core), also has little im- 
pact on wind anisotropies. Hydrodynamic simulations of the 
proto-neutron star formation and of the ac cretion during the 
first hundred s of m illiseconds after bounce dBuras et all |2006; 
iDessart et all 120061) show that the convective shell is placed at 
a density well above that of the neutrinosphere. The convec- 
tive layer is found to be separated from the neutrino spheric re- 
gion by a stable shell. The proto-neutron star convection in the 
post-explosion phase (later than one second after bounce) has 
not been explored by hydrodynamic simulations yet. We sus- 
pect — and conclude from the structure of our present models 
near the inner grid boundary — that this separation between con- 
vective layer and neutrinosphere is maintained also during the 
later stages. Our suspicion is justified by the decreasing optical 
depth in the vicinity of the neutrinosphere, which allows neu- 
trino transport to become more efficient than energy transport by 
convection. The convectively stable shell around the convective 
core of the proto-neutron star precludes convective asymmetries 
from having an immediate and undamped influence on the driv- 
ing layer of the wind outflow. Perturbations of the wind induced 
by convection thus remain small. 

iBuras et al.l (120061) : IDessart et aP (120061) showed that proto- 
neutron star convection speeds up the lepton number loss (by 
temporarily decreasing the electron- antineutrino luminosity) and 
increases the energy transport by higher muon and tau neutrino 
luminosities. At the same time it leads to reduced mean energies 
of the radiated neutrinos. While such effects are relevant for de- 
tailed modeling of the Kelvin-Helmholtz cooling of the nascent 
neutron star, they are subsumed by our treatment of the inner 
grid boundary. We therefore do not expect them to introduce any 
qualitatively new features into our models beyond what we can 
and have explored by variations of the boundary condition. 

Neutrino-flavor conversions near and above the proto- 
neutron star surface, i n particular due to neutrino-neutrino for- 
ward scattering (see iDuan et all I2Q1QL for a recent review) 
are disregarded in the present work, but may change the lu- 



minosity and spectral evolution of the different kinds of neu- 
trinos radiated by the forming neutron star. During the post- 
bounce accretion phase and the onset of the explosion, when 
the luminosities and spectra of different neutrino flavors show 
pronounced differences, such collective neutrino oscillations 
might have some relevance for the dynamics and for setting 
the neutron-to-proton ratio and thus the nucleosynthesis con- 
ditions in the early supernova ejecta. However, it is more un- 
certain whether any important effects remain during the Kelvin- 
Helmholtz cooling of the proto-neutron star after the explosion 
has been launched. Hydrodynamical calculations of this phase in 
spherical symmetry, using sophisticated neutrino transport meth- 
ods (but also ig n oring t he effects of n e utrino flavor oscillations; 
iHudepohl etaLl (l201Qh : iFischer et aP (120101) ) predict that after 
the end of the post-bounce accretion the luminosities and spectra 
for all kinds of neutrinos become very similar. Given that, flavor 
conversions will certainly have little influence on the neutrino- 
energy deposition and therefore on the dynamics of the neutrino- 
driven wind. Moreover, also their effect on the electron frac- 
tion in the neutrino-heated outflow is probably relatively minor. 
Nevertheless, subtle changes of the neutrino flux spectra due to 
energy- and direction-dependent oscillation effects cannot be ex- 
cluded. Since the neutron-to-proton ratio and the element forma- 
tion in the wind outflow are highly sensitive to details of the elec- 
tron neutrino and antineutrino luminosities and spectra, a closer 
investigation of this aspect may be desirable, but is very ambi- 
tious in the multi-dimensional context. 

3.4.3. General relativity 

The approximation of general relativistic gravity used in our 
simulations has been tested in great detail for its accuracy in 
different phases of stellar core collapse. C orresponding results 
were published in several preceding papers. iMarek et all (120061) 
compared the neutron star structure, stability properties, and os- 
cillation behavior obtained with the effective relativistic poten- 
tial and fully relativistic calcul ations of nonrotating and rotat- 
ing collapsing stellar iron cores. iLiebendor fer et al.l (120051) com- 
pared the dynamics and neutrino emission during the collapse 
phas e, core bounce, and early post-bounce evo lution in ID (see 
also IMarek et all 120061) . lArcones et aP (120071) (Paper I) inves- 
tigated the structure of the neutrino-driven win d in compari- 
son to steady- state rel ativistic wind solut ions. iHiidepohl et al.l 
( 2010) and in particular Mu ller et aP (1201 Ol) explored the Kelvin- 
Helmholtz cooling evolution of proto-neutron stars and their 
neutrino emission. In these studies fully relativistic models were 
contrasted with simulations in which approximated relativis- 
tic gravity (plus gravitational redshift in the neutrino transport 
where applied) were combined with otherwise Newtonian hy- 
drodynamics. 

The core-collapse phase, core bounce, shock break-out, and 
early accretion were found to be in excellent agreement, the 
neutron-star structure, its binding energy, neutrino emission 
properties and total energy loss were found to be reproduced in 
all quantities to better than 10% and in most cases and phases 
even to percent accuracy, and the neutrino-driven wind profiles 
turned out to be nicely compatible with steady- state relativistic 
solutions. Although there can be differences in smaller details, 
the main effects of general relativity seem to be captured well by 
our approximative treatment. The reasons for this result are the 
facts that (1) the effective relativistic potential was constructed 
to explicitly include the relativistic modifications of the Tolman- 
Oppenheimer-Volkoff eq uations compared to Newtonian hydro- 
static equilibrium (see lMarek et all 120061) . and (2) the fluid ve- 
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locities in the supernova core hardly ever exceed 10-20% of the 
speed of light, in which case Newtonian fluid dynamics is still a 
good approximation. 

3.4.4. Neutron star motion 

A kicked neutron star will move away from the grid center but 
will remain the center of the neutrino-driven wind. The neutron 
star motion may thus lead to a distortion of the wind-ejecta in- 
teraction and may change the location of the reverse shock in the 
wind and thus its influence on the wind properties. Such influ- 
ence on the detailed e volution of a particular model was indeed 
seen by lScheck et al.l (120061) . who performed runs (for fixed pro- 
genitor, boundary conditions, and perturbation seeds) without 
and with applying a Galilei transformation on the flow around 
a neutron star at the grid center by giving the surrounding fluid 
a coherent motion with the negative sign of the neutron star kick 
velocity. This procedure was assumed to capture the main effects 
of the neutron star movement. 

Although such effects can be relevant for determining the 
detailed evolution and ejecta properties of a particular model star 
(which was not the goal of this work) they are unlikely to change 
the basic and general conclusions that we have drawn from our 
simulations. Since the neutron- star velocity v ns is typically small 
compared to the sound speed and the terminal wind velocity, 
Vns ^ £s and v ns <*c u w (r » R m ), it is hard to imagine that 
allowing the neutron star to move (instead of anchoring it at the 
grid center as in our models) will affect the overall picture that 
we have obtained for the anisotropic wind-ejecta interaction. 

3.4.5. Three-dimensional modeling 

Performing simulations in two dimensions imposes the artifical 
constraint that all structures are axially symmetric. Of course, 
this naturally raises the question what differences one might ex- 
pect in three dimensions. 

Ignoring the important question whether 2D/3D differences 
have any influence on the success of the neutrino-heating mech- 
anism for supernova explosions ( a very first statement i n this 
context was recently published by Nordha us et al.l (120101) ). and 
relying on the viability of this mechanism (which is a funda- 
mental assumption in our study), basic features of the 2D asym- 
metries observed in ou r calculations seem to be co nfirmed by 
recent 3D simulations: IWongwathanarat etaD (1201 Oh followed 
the evolution of a small set of 3D models over a similarly long 
post-bounce period (but with considerably less spatial resolu- 
tion than used in our 2D calculations). The 3D slice of figure 3, 
right panel, in the latter publication exhibits the same overall be- 
havior and structural features for the neutrino-driven wind phase 
in 2D and 3D models, namely an essentially spherical wind, an 
asymmetric wind-termination shock, and an inhomogeneous and 
strongly anisotropic ejecta shell behind the forward shock. Also 
in 3D the reverse shock is highly deformed and its effect on the 
wind properties and nucleosynthesis conditions depends on the 
outflow direction. We are therefore confident that our main find- 
ings for the wind-ejecta interaction are valid not only in the con- 
sidered 2D situation. 

Of course, the present parameterized explosion models, 
which do not yield any information about the explosion prop- 
erties (energy, timescale) of a progenitor star, cannot be conclu- 
sive on the consequences of these asymmetries for the explosive 
nucleosynthesis of individual stars, nor does our present knowl- 



edge of the explosion mechanism allow for any statements in a 
statistical sense. 

Rotation of the nascent neutron star is an additional degree 
of freedom, which we ignored in the models of this paper. It will 
cause a global pole-equator asymmetry of the wind and of the 
wind-ejecta interaction with potentially interesting implications. 
This should be studied in future work by systematic variations 
of the proto-neutron star spin. 

4. Nucleosynthesis implications 

We have shown the impact of multidimensional effects on the 
dynamical evolution of the neutrino-driven wind, reverse shock 
and supernova ejecta. In this section, we want to briefly ad- 
dress the possible implications of our results for the nucle- 
osynthesis processes occurring in s upernova outflows: charged- 
particle reactions , alpha process (Wooslev & Hoffman, 1992; 
Witt i et all 11994. vp -process dFrohlich et all 120061: iPruet et all 
2006; Wanajo, 2006), and occasionally r-process (lArnould et all 
|2007[ for a review). 

Sin ce the works of ICameronl (Il957l) and Burbidg e et al.l 
(fl957h . core-collapse supernova outflows have been the best 
studied candidate for the production of heavy elements. 
However, this environment is facing more and more problems 
to fulfill the requirements (high entropy, low electron fraction 
and fast expansion) for the production of heavy r-process el- 
ements (A> 90). The conditions found to be necessa ry for a 
robust and strong r-process (e.g., Hoff manet all 1 19971) are not 
achieved by recent long-time supernova simulations (Paper I, 
Hiid epohl etalll2010l:lFischer et al.Ll201Qb . This is also the case 
for our 2D simulations, where the wind entropies are too low to 
get the high neutron-to-seed ratios necessary for the r-process. 

Yet galac ti c chemical evolution m od els ( see e.g., 
Ilshimaruetail 12004 lOian & WasserburgL 120081) suggest 
that at least a subset of core-collapse supernovae could be 
responsible of the origin of half of the heavy r-process elements. 
Therefore, one may speculate that the r-process could take 
place in neutrino-driven winds because of some still unknown 
aspect of physics that might cure the problems revealed by 
the present hydrodynamical models. In this case the revers e 
shock could have important consequences (Wan ajo et al.Ll2002b . 
Depending on the temperature at the reverse shock the r-process 
path is different. When the reverse- shock temperature is low 
(r rs < 0.5 GK), neutron-capture an d beta-decay timescales 
are similar (iBlake & Schramm! Il976l) and shorter tha n (y, n) 
timescales. This is also kn own as "cold r-process" (Wanaio, 
2007; Panov & Janka, 2009). In contrast, when the reverse shock 
is at high temperatures, t here is an (n , y) - ( y, n) equilibrium as 
in the classical r-process dKratz et all 1 19931: iFreiburghaus et all 
1 19991: iFarouqi et all l2010h . The final abundances for these two 
types of evolution are very different , see e.g .. IWanaio et al. 
(2002); Wanaio (2007 ^ lKurodaetal.1 (120081) : iPanov & Janka 
(200§; lArcones & Martinez-Pinedol (l2010h . Also the impact 
of the n uclear physics input varies depend ing on how matter 
expands (lArcones & Martinez-Pinedol 1201 Oh . 

In addition to the r-process, whose astrophysical site is 
still uncertain, there are other nucleosynthesis processes occur- 
ring in supernova outflows. The reverse shock can affect the 
production of p-nuclei that happ ens in neutrino-driven w inds 
via charged-pa rticle reactions (IWooslev & Hoffman, 1992) and 
the v p-process (Frohli chet all 12006: Pru et et aLU 2006: Wanajo, 
2006). This is becoming more important, because the most re- 
cent and most sophisticated supernova simulati ons show that 
the ejecta are proton-rich for several seconds (iFischer et all 
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2010) and even until com pletion of the proto-n eutron star cool- 
ing and deleptonization (Hude pohl et all 120101) . The impact of 
the reverse shock on these nucleosynthesis processes can show 
up in two ways. First, the temperature jump leads to an in- 
crease of the photo-dissociation rate. When the latter is too high, 
newly formed nuclei are destroyed. However, when the photo- 
dissociation rate is moderate, the temperature increase favors 
the captures of charged particles. The other important effect of 
the reverse shock is the strong reduction of the expansion ve- 
locity with the consequence that the temperature stays constant 
or decreases only slowly after a mass element has crossed the 
reverse shock. Depending on the exact value of the tempera- 
ture, photo-dissociations or charged-particle reactions continue 
to take place. Moreover, when the expansion becomes slower, 
the matter stays exposed to high neutrino fluxes for a longer 
time. This increases the efficiency of the vp-process. However, 
one should notice that the processes described here are possi- 
ble only when the reverse- shock radius is sufficiently small, e.g., 
during the first few seconds after the onset of the explosion, be- 
cause otherwise the temperature at t he reverse sho c k is al ready 
too low to play any role. Recently, IWanajo et al.l (120101) sug- 
gested a possible significant impact of the reverse shock on the 
nucleosynthesis un der these conditions, in contrast to the small 
effects reported by lRoberts et al.l (l2010h . Therefore, further nu- 
cleosynthesis studies should be done, taking into account the 
reverse- shock behaviour found in Paper I and in particular in the 
present work, where a huge variability due to multi-dimensional 
effects was obtained. 

5. Conclusions 

With a small set of two-dimensional simulations for three pro- 
genitor stars of different masses (Table [B we have demon- 
strated that the reverse- shock radius and the conditions of the 
shocked neutrino-driven wind matter in supernova explosions 
become strongly angle dependent. As we found in Paper I by 
analytic means, the position of the reverse shock depends on 
the wind properties (mass outflow rate and velocity) as well as 
on the pressure of the more slowly moving supernova ejecta. 
Comparison of the radial profiles of ID and 2D simulations 
shows that the neutrino-driven wind is spherically symmetric. 
This is caused by the isotropic neutrino emission from a neu- 
tron star that stays spherical in the absence of rotation. Multi- 
dimensional simulat ions including rotation could lead to differ- 
ences in the wind dMetzger et al.U2007l:IWanaioU2006l) and thus 
significant changes in its interaction with the supernova ejecta. 
Although without rotation the wind develops identically in all di- 
rections, a strong angular variation of the reverse shock position 
can appear because of the anisotropic matter (and thus pressure) 
distribution in the more slowly moving, early supernova ejecta. 

When the radial location of the reverse shock is constrained 
by the existence of dense downdrafts in the ejecta shell that fol- 
lows the supernova shock, the angle between the reverse shock 
and the wind velocity, which goes in the radial direction, can 
be estimated. We have found an analytic expression that relates 
this angle to the jump of the pressure at the reverse shock. The 
presence of the downflow features in the supernova ejecta shell 
with local density and pressure maxima leads to angular vari- 
ations of the radius of the reverse shock and thus a deforma- 
tion of the wind- shell boundary with kinks appearing in the re- 
gions where the obliqueness of the shock abruptly changes. The 
Rankine-Hugoniot conditions imply that oblique shocks are less 
efficient in decelerating matter. This produces a collimation of 
the shocked flow in the vicinity of the kinks. Finally, we have 



proven that basic features of the progenitor dependence seen in 
Paper I are also present in our two-dimensional simulations, e.g., 
a slower expansion of the forward and reverse shocks in more 
massive stars. However, such general aspects are superimposed 
by an enormous amount of variability of explosions even of the 
same progenitor and similar explosion energy due to the chaotic 
growth of nonradial hydrodynamic instabilities from small ini- 
tial seed perturbations. 

In summary, we have found that in the multi-dimensional 
case the expansion of the wind matter varies with the angular di- 
rection because it is influenced by the anisotropic distribution of 
the earlier ejecta, which evolves chaotically from initial random 
perturbations. Therefore, we strongly recommend that future nu- 
cleosynthesis studies should test the effect of different extrapo- 
lations of the evolution of the shock-decelerated ejecta. Our re- 
sults can be used as guidance for the overall variability that is 
possible and affects the nucleosynthesis-relevant conditions in 
multi-dimensional supernova environments. 
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Fig. 6. Profiles of radial velocity, pressure, density, temperature, 
and entropy as functions of radius at time 1.5 s after bounce. The 
one-dimensional model M10-ll-rl (black lines) is compared to 
profiles of the two-dimensional model T10-ll-rl at angles 9 = 
25 degrees (green lines) and 6 = 100 degrees (red lines). 
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Fig. 7. Evolution of reverse shock radius, and of the entropy, 
pressure, density, and temperature of the shocked neutrino- wind 
matter for models M10-ll-rl (black lines) and T10-ll-rl at = 
25 degrees (green lines) and 6 = 100 degrees (red lines). 



